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The Age of the Universe 
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Globular clusters are the oldest objects in the Galaxy whose age may be accurately determined. As such globular 
cluster ages provide the best estimate for the age of the universe. The age of a globular cluster is determined 
by a comparison between theoretical stellar evolution models and observational data. Current uncertainties in 
the stellar models and age dating process are discussed in detail. The best estimate for the absolute age of the 
globular clusters is 14.6±1.7 Gyr. The one-sided, 95% confidence limit on the lower age of the universe is 12.2 Gyr. 
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1. Introduction 

A minimum age for the universe may be deter- 
mined by obtaining a reliable estimate for the age 
of the oldest objects within the universe. Thus, in 
order to estimate the age of the universe (t Q ), the 
oldest objects must be identified and dated. The 
oldest objects in the universe should contain very 
little (if any) heavy elements as nucleosynthesis 
during the Big Bang only produced hydrogen, 
helium, and lithium. All elements heavier than 
lithium are not primordial in origin, and their 
presence indicate that an object was not the first 
to have formed in the universe. Unfortunately, 
astronomers have been unable to locate any ob- 
ject which consists solely of primordial elements. 
There are, however, objects which contain small 
amounts of the heavier elements. These objects 
will be the focus of my talk. 

Within our galaxy, the oldest objects which can 
be dated are the globular clusters (GCs). GCs are 
compact stellar systems, containing ~ 10 5 stars 
(see figure 1). These stars contain few heavy ele- 
ments (typically 1/10 to 1/100 the ratio found in 
the Sun). GCs are spherically distributed about 
the Galactic center, suggesting that GCs were 
formed soon after the proto-Galactic gas started 
collapsing. Thus, GCs were among the first ob- 
jects formed in the Galaxy. An estimate of their 
age will provide a reasonable lower limit to the 
age of the universe. In order to estimate the true 
age of the universe, one must add to the age of 
the GCs the time it took GCs to form after the 
Big Bang. 




Figure 1. An example of a typical GC, M53. This 
cluster contains roughly 10 6 stars. Image cour- 
tesy of Ata Sarajedini. 

It is important to realize that the estimate for 
t a is based on a single method, applying stellar 
evolution models to GC observations. This is in 
sharp contrast to estimates for the Hubble con- 
stant (Ho), which are based on a wide variety of 
independent techniques (see the review by Trim- 
ble in this volume). Whereas estimates for the 
Hubble constant may differ by a factor of two, de- 
pending on the observer and the technique which 
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is used, estimates for the absolute age of the GCs 
typically agree with each to within ~ 10%. The 
important considerations then, are to estimate 
the uncertainty in the stellar models and the de- 
rived value of t , and to test the stellar models 
in as many ways as possible to ensure that no 
systematic errors exist. 

The study of t a has taken on increasing signif- 
icance in recent years, as the value of t a derived 
from GCs appears to be longer than the expan- 
sion age of the universe, derived from H . The 
value of H t is a function of the the cosmologi- 
cal constant (A) and density of the universe (f2, 
in units of the critical density). For a 'standard' 
inflationary universe (fl = 1, A = 0), H t = 2/3. 
A value of H t > 1 requires a non-zero cosmolog- 
ical constant, or a significant revision to standard 
Big Bang cosmology. 

In this review I will describe in detail how GC 
ages are estimated. Some basic observational 
properties of GCs are summarized in §|^. The con- 
struction of stellar models which are used to date 
GCs are described in §[| while §jj contains a dis- 
cussion of age determination techniques for GCs. 
Section contains a detailed discussion of possi- 
ble errors in the age estimates for GCs. This sec- 
tion includes the results of a recent Monte Carlo 
analysis which has resulted in a firm lower limit 
to the age of the universe. Various tests of stel- 
lar models are presented in including white 
dwarf cooling time-scales. Finally, §0 contains a 
summary of this review. 

2. Globular Cluster Observations 

Observers typically measure the apparent mag- 
nitude (oc 2.5 \og(luminosityj) of as many stars 
as possible within a GC. These measurements are 
usually taken through at least two different fil- 
ters, so that the apparent colour of the stars may 
also be determined. The fact that a GC contains 
a large number of stars all at the same distance 
from the Earth is an enormous advantage in in- 
terpreting the observations. The ranking in ap- 
parent luminosity (as seen in the sky) is identical 
to the ranking by absolute luminosity. Figure 2 
is a typical example of how the observations are 
reported, as a colour-magnitude diagram. 
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Figure 2. A colour- magnitude diagram of a typical 
GC, M15 [1]. The vertical axis plots the magni- 
tude (luminosity) of the star in the V wavelength 
region, with brighter stars having smaller magni- 
tudes. The horizontal axis plots the colour (sur- 
face temperature) of the stars, with cooler stars 
towards the right. The various evolutionary se- 
quence have been labeled (see text). For clarity, 
only about 10% of the stars on the main sequence 
have been plotted. 

One of the great triumphs of the theory of stel- 
lar structure and evolution has been an explana- 
tion of the colour-magnitude diagram. The more 
massive a star, the quicker it burns its nuclear fuel 
and evolves. Thus, stars of different initial masses 
will be in different stages of evolution. Figure 2 
graphically illustrates the various phases of stellar 
evolution. After their birth, stars start to burn 
hydrogen in their core. This is referred to as the 
main sequence. A star will spend approximately 
90% of its life on the main sequence. The Sun is 
a typical example of a main sequence star. Even- 
tually, a star exhausts the supply of hydrogen in 
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its core, and reaches the main sequence turn-off 
(MSTO). This point is critical in the age deter- 
mination process. After a star has burned all of 
the hydrogen in its core, the outer layers expand 
and hydrogen fusion occurs in a shell surround- 
ing the helium core. The expansion of the outer 
layers causes the star to cool and become red, so 
stars in this phase of evolution are said to occupy 
the red giant branch. The hydrogen burning shell 
moves out in mass coordinates, leading to increas- 
ing luminosity and helium core mass. On the red 
giant branch, a typical GC star is believed to lose 
ss 25% of its mass. When and how this occurs 
is a still a subject of research. Eventually, the 
helium core becomes so dense that helium fusion 
is ignited. The star quickly settles onto the hori- 
zontal branch (HB). On the HB, fusion of helium 
occurs in the core, surrounded by a shell of hydro- 
gen fusion. Exactly where a star lies on the hori- 
zontal branch (blue or red) depends on how much 
mass loss has occurred on the red giant branch. 
Some stars on the horizontal branch are unsta- 
ble to radial pulsations — these stars are referred 
to as RR Lyrae stars. A star's lifetime on the 
HB is extremely short, it soon exhausts the sup- 
ply of helium at its core and becomes an asymp- 
totic giant branch star (similar to the red giant 
branch), burning helium and hydrogen in shells 
about a carbon core. In GC stars (like the Sun), 
the core temperatures and densities never become 
high enough to ignite the fusion of carbon. Af- 
ter a star finishes its helium and hydrogen shell 
burning, the envelope may be ejected, while the 
core contracts and becomes extremely dense. The 
star becomes dim, as the only energy available to 
the star is that from gravitational contraction. 
In this terminal phase of evolution, the star is re- 
ferred to as a white dwarf (not shown in Fig. 2). 
When the Sun becomes a white dwarf, its radius 
will be similar to the Earth's. As a white dwarf 
ages, it continues to cool, and emit less radiation. 
Ultimately, a star will reach equilibrium with its 
surroundings, becoming virtually invisible. 

3. Stellar Models 

Our understanding of stellar evolution is based 
on stellar structure theory. There are numerous 



textbooks which describe the basic theory of stel- 
lar structure and the construction of stellar mod- 
els (e.g. 0). A stellar model is constructed by 
solving the four basic equations of stellar struc- 
ture: (1) conservation of mass; (2) conservation 
of energy; (3) hydrostatic equilibrium and (4) en- 
ergy transport via radiation, convection and/or 
conduction. These four, coupled differential equa- 
tions represent a two point boundary value prob- 
lem. Two of the boundary conditions are speci- 
fied at the center of the star (mass and luminosity 
are zero), and two at the surface. In order to solve 
these equations, supplementary information is re- 
quired. The surface boundary conditions (tem- 
perature and pressure) are based on stellar atmo- 
sphere calculations. The equation of state, opac- 
ities and nuclear reaction rates must be known. 
The mass and initial composition of the star need 
to be specified. Finally, as convection can be im- 
portant in a star, one must have a theory of con- 
vection which determines when a region of a star 
is unstable to convective motions, and if so, the 
efficiency of the resulting heat transport. Once all 
of the above information has been determined a 
stellar model may be constructed by solving the 
four stellar structure equations. The evolution 
of a star may be followed by computing a static 
stellar structure model, updating the composition 
profile to reflect the changes due to nuclear reac- 
tions and/or mixing due to convection, and then 
re-computing the stellar structure model. 

There are a number of uncertainties associated 
with stellar evolution models, and hence, age es- 
timates based on the models. Probably the least 
understood aspect of stellar modeling is the treat- 
ment of convection. The understanding of convec- 
tion in a compressible plasma, where significant 
amounts of energy can be carried by radiation, is 
a long standing problem. Numerical simulations 
hold promise for the future Q , but at present one 
must view properties of stellar models which de- 
pend on the treatment of convection to be uncer- 
tain, and subject to possibility large systematic 
errors. Main sequence, and red giant branch GC 
stars have surface convection zones. Hence, the 
surface properties of the stellar models (such as 
its effective temperature, or colour) are rather un- 
certain. Horizontal branch stars have convective 
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cores, so the predicted luminosities and lifetimes 
of these stars are subject to possible systematic 
errors. 

Another important consideration in assessing 
the reliability of stellar models, and the ages they 
predict for GCs is that the advanced evolutionary 
stages are considerably more complicated than 
the main sequence. Thus, one may expect that 
the main sequence models are least likely to be 
in error. Observations of CNO abundances in 
red giant branch stars indicate that some form 
of deep mixing occurs in these stars, which is not 
present in the models Q. In contrast, there is 
no observational evidence suggesting that the low 
mass, main sequence models are in serious error. 
For this reason, age indicators which are based on 
main sequence models are the most reliable. 

4. Globular Cluster Age Estimates 

A theoretical stellar evolution model follows 
the time evolution of a star of a given initial 
composition and mass. Stars in a given GC 
all have the same chemical composition and age 
(with a few exceptions), but different masses. 
Thus, in order to determine the age of a GC, 
a series of stellar evolution models with the 
same composition but different masses must be 
constructed. Interpolation among these mod- 
els yields an isochrone, a theoretical locus of 
points for stars with different masses but the same 
age. The theoretical calculations are performed 
in terms of total luminosity, and effective tem- 
perature. As discusses in observers measure 
the brightness of a star over specified wavelength 
ranges. Thus, it is necessary to convert from the 
effective temperature and total luminosity to lu- 
minosities in a few specified wavelength intervals. 
This requires a detailed knowledge of the pre- 
dicted flux, as a function of wavelength. Theoreti- 
cal stellar atmosphere models are used to perform 
this conversion. The result of such calculations is 
illustrated in Figure 3, which plots isochrones for 
3 different ages, in terms of absolute V magni- 
tude, and B-V colour. As can be seen in Figure 
3, differences in age lead to large differences in 
the MSTO region. The MSTO becomes fainter, 
and redder as the age increases. 
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Figure 3. Theoretical isochrones for 3 ages (10, 14 
& 18 Gyr). Note that the main sequence turn-off 
becomes fainter and redder for older isochrones. 



In order to provide the tightest possible con- 
straints on the age of the universe, it is important 
to use an age indicator which has the smallest 
possible theoretical error. From the discussion in 
§[| it is clear that (a) the main sequence is the 
best understood phase of stellar evolution and (b) 
the predicted luminosities of the models are bet- 
ter known than the colours. These two reasons, 
coupled with the age variation exhibited in Figure 
3 lead to the conclusion that the absolute mag- 
nitude of the MSTO is the best indicator of the 
absolute age of GCs || . Unfortunately determin- 
ing the magnitude of the MSTO in observational 
data is quite difficult, as the MSTO region ex- 
tends over a large range in magnitude (see Fig. 
2). For this reason, it is best to determine the 
mean age of a large number of GCs, in order to 
minimize the observational error. 
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Observers measure the apparent magnitude of 
a star. In order to convert to the absolute mag- 
nitude (and hence, determine an age), the dis- 
tance to a GC must be determined. Obtaining 
the distance to an object remains one of the most 
difficult aspects of astronomy. At present, there 
are two main techniques which are used to deter- 
mine the distance to a GC (I) main sequence fit- 
ting to local sub-dwarfs with well measured par- 
allaxes, and (2) using the observed magnitude of 
the HB combined with a relationship for the ab- 
solute magnitude of the HB (derived using RR 
Lyrae stars). Unfortunately, there are few sub- 
dwarfs with well measured parallaxes (a situation 
which should change once date from the Hippar- 
cos satellite are released), so at the present time 
the use of the HB to set the distance scale to 
GCs is the most reliable. The HB has the ad- 
vantage that the difference in magnitude between 
the MSTO and the HB (AV(TO - HB)) is inde- 
pendent of reddening. Thus, AV(TO — HB) is a 
widely used age determination technique, which 
uses the absolute magnitude of the main sequence 
turn-off as its age diagnostic. There are a number 
of observational and theoretical techniques which 
may be used to obtain the absolute magnitude 
of the RR Lyr stars (M V (RR)), with the general 
conclusion that M V (RR) = /i [Fe/H] + 7 where /z 
is the slope with metallicity and 7 is the zero- 
point Q . Uncertainties in the slope primarily ef- 
fect the relative ages of GCs. When a number of 
GCs are studied, the uncertainty in the slope has 
a negligible effect on the derived mean age. The 
uncertainty in the M v (RR) zero-point has a large 
impact on the derived ages (see §||). 

A number of different researchers have 
constructed stellar evolutionary models and 
isochrones which they have used to estimate the 
age of the GCs. These estimates agree well with 
each other, and indeed have remained relatively 
constant since ~ 1970 This is not too surpris- 
ing, as the basic assumptions and physics used to 
construct the stellar models are the same for the 
different research groups, and have not changed 
for a number of years. These studies have also re- 
vealed that GCs are not all the same age. Thus, 
to provide the best estimate for the age of the 
universe one must select a sample of the oldest 



GCs. Chaboyer, Demarque, Kernan & Krauss 
(hereafter CDKK, ||]) have recently completed a 
study of the absolute age of the oldest GCs, which 
they found to be 14.6 Gyr. The novel aspect of 
this work was the detailed consideration of the 
possible sources of error in the stellar models and 
age determination process which allowed CDKK 
to provide an estimate of the error associated with 
their age estimate. 

5. Error Estimates 

To access the error in the absolute GC age es- 
timates, one must review the assumptions and 
physics which are used to construct stellar mod- 
els and isochrones. The discussion presented in 
§H and §|| of the GC ages determination process 
allows one to compile a list of possible sources of 
error in theoretical calibration of age as a function 
of the absolute magnitude of the MSTO: 

1. assumption of hydrostatic equilibrium in ra- 
diative regions 

2. nuclear reaction rates 

3. opacities 

4. equation of state 

5. treatment of convection 

6. surface boundary conditions 

7. chemical composition 

8. conversion from theoretical luminosities to 
observed magnitudes 

This is a lengthly list, which has been studied in 
some detail In this review, I will concen- 

trate on a few items which turn out to be partic- 
ularly important, or for which improved calcula- 
tions have recently become available. 

The validity of the assumption of hydrostatic 
equilibrium in the radiative regions of stars has 
received considerable attention. Not surprisingly, 
if some process operates which mixes material 
into the core of a star, the main sequence life- 
times will be prolonged, and hence the true age 
of the GCs will be older than current estimates. 



G 



However, if a microscopic diffusion is active (caus- 
ing helium to sink relative to hydrogen) , than the 
main sequence life-times are shortened, leading 
to lower estimates for the age of the GCs |Io| . 
The inclusion of diffusion has been found to lower 
the GC age estimates by 7%. There is evidence 
from helioseismology that diffusion is occurring 
in the sun but models of halo stars which 
incorporate diffusion are unable to match the Li 
observations in halo stars |l2|. Until this mat- 
ter is resolved, GC age estimates are subject to a 
possible 7% systematic error due to the effects of 
diffusion. 

The correctness of the nuclear reaction rates 
used in stellar models have been extensively an- 
alyzed due to their importance in the solar neu- 
trino problem jl3| . Although predicted solar neu- 
trino fluxes are quite sensitive to possible errors 
in the nuclear reaction rates, they have a minor 
effect on GC age estimates |J . 

There has been considerable effort devoted to 
determining the opacities used in stellar mod- 
els. A number of different research groups have 
calculated opacities, using independent methods. 
The agreement between these calculations, par- 
ticularly for metal-poor mixtures is quite good 
Indeed, the high-temperature opacities in 
metal-poor stars have only changed by ~ 1% be- 
tween the mid-1970's and the present. 

The equation of state used in stellar models 
is another area of active research. Detailed cal- 
culations have lead to the availability of equa- 
tion of state tables which are a considerable im- 
provement over the simple analytical formulae 
usually employed in stellar evolution calculations 
pq |. These calculations include Coloumb effects, 
which had typically been ignored in previous cal- 
culations. It has been found that the improved 
equation of state reduces GC age estimates by 7% 
p6[ | . Independent calculations (using entirely dif- 
ferent physical assumptions) |17| lead to an equa- 
tion of state which agrees quite well with |15| . 
Thus, it is unlikely that there are significant er- 
rors in the new equation of state calculations. 

The correct composition to use in stellar mod- 
els is a long standing problem. The helium mass 
fraction is taken to be the primordial value Jl8| , 
Y = 0.23 — 0.24. A generous ±0.03 uncertainty in 



Y leads to a negligible uncertainty in the derived 
age ||. In contrast, the uncertainty in the heavy 
element composition leads to a significant error 
in the derived age. It is relatively easy to deter- 
mine the abundance of iron, and it was generally 
assumed that the other heavy elements present 
in a GC star would be scale in a similar manner. 
However, it has become clear that the a-capture 
elements (O, Mg, Si, S, and Ca) are enhanced in 
GC stars, relative to the ratio found in the Sun 
fT9[] . Oxygen is the most important of these ele- 
ments, being by far the most abundant. However, 
it is quite difficult to determine the abundance of 
oxygen in GC stars (20) , and as a consequence the 
uncertainty in the oxygen abundance is leads to 
a relatively large error (±6%) in the derived age 
of the oldest GCs @. 

In an attempt to take into account all of the 
possible uncertainties in the GC age dating pro- 
cess, CDKK have made a detailed examination of 
the likely error in each of the parameters. CDKK 
performed a Monte-Carlo simulation, in which 
the various quantities were allowed to vary within 
some specified distribution, chosen to encompass 
the possible uncertainty in that quantity. The 
mean age of the 17 oldest GCs was determined 
using 1000 independent sets of isochrones. As- 
suming that the distances to the GCs are known 
exactly, a mean age of 14.6 ± 1.1 Gyr was deter- 
mined, with a 95% confidence limit that GCs are 
older than 12.9 Gyr. Allowing for an uncertainty 
in the distance modulus (±0.08 mag in the RR 
Lyrae zero-point) increases the allowed range to 
14.6 ± 1.6 Gyr, with a 95% confidence limit that 
the GCs are older than 12.1 Gyr. The error in 
the distance modulus dominates the overall un- 
certainty in the absolute age of the oldest GCs. 

6. Tests of Stellar Models 

The discussion presented in the previous sec- 
tion, and the error analysis performed by CDKK 
assumes that there are no unknown systematic 
errors in the GC age determination process. This 
assumption is supported by four independent 
tests of stellar structure theory: 

1. fitting theoretical isochrones to observed 
GC colour-magnitude diagrams; 
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2. comparison between observed and predicted 
luminosity functions; 

3. observations of solar p-modes which probe 
the structure of the Sun down to r = 
0.05 R ; and 

4. white dwarf age estimates for the local 
Galactic disk which agree with MSTO age 
estimates. 

These tests are summarized below. 

Theoretical isochrones provide a good match 
to observed GC colour-magnitude diagrams (Fig- 
ure 4). The freedom to modify the predicted 
colours of the models (due to our poor treatment 
of the surface convection in these stars), implies 
that this is not a definitive test of stellar mod- 
els which proves that they are correct. However, 
the absence of any unexplained features in the ob- 
served colour-magnitude diagram constrains non- 
standard models. For example, models which in- 
clude a mixed core (which predict older ages for 
GCs) predict a 'hook' in the MSTO region which 
is not observed. Hence, one may conclude that 
GC stars do not have cores which have been ex- 
tensively mixed. 

The number of stars as a function of luminos- 
ity is referred to as a luminosity function (LF). 
On the lower main sequence, the LF is a func- 
tion of the number of stars per unit mass, and 
the mass-luminosity relationship. In the more ad- 
vanced evolutionary stages (starting about 1 mag- 
nitude below the MSTO), the evolutionary time- 
scales are very short, and dominate the number 
counts. Hence, observed LFs provide an excel- 
lent test of the relative lifetimes predicted by the 
stellar models. The freedom to choose an over- 
all normalization factor between the observations 
and theory implies that this is not a test of the 
absolute lifetimes. In general, a good match is 
found between predicted and observed LFs |^2) 
(see Figure 5), implying that the relative evolu- 
tionary time-scales predicted by the models are 
correct. Thus, any mechanism which shortens the 
main sequence lifetime of GC stars (and hence, 
shortens the GC age estimates) must predict a 
corresponding decrease in the more advanced evo- 
lutionary phases, like the RGB. There are sugges- 
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Figure 4. Theoretical isochrones with ages of 14, 
16 and 18 Gyr are fit to the GC M68 [21]. These 
isochrones do a good job of matching the observed 
colour-magnitude diagram. 



tions with the present data sets that models which 
incorporate isothermal cores do not match the ob- 
servations p^ ], although the conclusions are not 
definitive. The isothermal core models predict 
GC ages which are about 20% smaller than the 
standard stellar models [Q . It is now technically 
possible to obtain much larger observational data 
sets, which will lead to much smaller (Poisson) er- 
ror bars, and a more definitive test of the relative 
lifetimes predicted by stellar evolution models. 

The internal structure of the Sun is predicated 
to be quite similar to a typical main sequence 
GC star. Both stars fuse hydrogen via the pp 
cycle in their radiative interiors and have a sur- 
face convection zone. Hence, tests of solar models 
may also be viewed as tests of the stellar models 
which are used to determine GC ages. Millions of 
non-radial oscillatory modes have been observed 
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Figure 5. Theoretical luminosity functions with 
ages of 12, 14 and 16 Gyr are fit to the GC NGC 
288 [23]. The luminosity functions do a good job 
of matching the observed data. 



at the surface of the Sun. These non-radial modes 
are referred to as p-modes, and provide an unique 
test of stellar evolution. Precise observations of 
the frequencies of the p-modes make it possible 
to infer many properties of the solar interior and 
to test stellar evolution models p5| . For example, 
helioseismology has lead to estimates for the solar 
helium abundance |26| and has put strict limits on 
the amount of overshoot present below the surface 
convection zone p7| . Inversions of the observed 
frequencies with respect to a solar model yield the 
difference in the squared sound speed between the 
model and the Sun. The squared sound speed is 
proportional to the pressure divided by the den- 
sity, c 2 oc P/p. Thus, inversions of the solar p- 
modes are a direct test of the interior structure 
predicted by solar models. The results of such an 
inversion are shown in Figure 6. The agreement 
is remarkable good, with differences of less than 
0.5% throughout most of the model. The inver- 
sions do not extend to the very center of the star 
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Figure 6. The difference in the square of the 
sound speed, between a reference solar model and 
the actual Sun, as a function of radius. This plot 
is obtained by inverting the observed frequencies 
of the solar p-modes with respect to the reference 
model. 



- the observed p-modes do not penetrate below 
r ~ 0.05 R Q , and so do not probe the structure 
of the sun below this point. The p-mode obser- 
vations indicate that the surface structure of the 
models are in error (implying that the treatment 
of convection needs to be improved). 

The excellent agreement between the sound 
speed in the Sun, and that predicted by solar 
models is strong evidence that there are no se- 
rious errors in current stellar evolution models. 
However, there remains the long standing discrep- 
ancy between the predicted solar neutrino fluxes 
and those observed on the earth. Four indepen- 
dent neutrino experiments have observed a solar 
neutrino flux which is 1/2 — 1/3 the predicted 
value [^8). A solution to this problems requires 
either (a) new neutrino physics, or (b) a system- 
atic error in the stellar evolution models. Given 
the excellent agreement with helioseismology, and 
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the apparent energy dependence of the observed 
solar neutrino deficient, it is likely that a resolu- 
tion of the solar neutrino problem requires new 
neutrino physics |29[) . However, until definitive 
observational evidence is obtained Q, there re- 
mains a possibility that there is some unknown 
systematic error in the solar models. If this were 
to be the case, then our estimates for GC ages 
would require revision. 

Estimates for the age of a stellar popula- 
tion may be obtained from white dwarf cooling 
curves. The assumptions and physics used to 
construct white dwarf cooling models are quite 
distinct from those used in stellar evolution mod- 
els. Hence, white dwarf cooling curves provide 
an independent test of the lifetimes predicted by 
stellar evolution models. The basic idea behind 
white dwarf cooling curve age estimates is that 
as a white dwarf ages, it becomes fainter, and 
cooler. Thus, for a given age white dwarfs will 
not exist below a minimum temperature and lu- 
minosity. White dwarf cooling curves are rela- 
tively simple to model, but it is difficult to observe 
white dwarfs due to their low luminosity. Indeed, 
at the present time it is impossible to detect the 
faint end of the white dwarf luminosity function 
in GCs. The turn-over in the white dwarf lu- 
minosity function has been detected in the local 
solar neighborhood, and provides an independent 
estimate for the lifetime of the Galactic disk of 
10. Gyr This is in agreement with esti- 
mates for the age of the oldest open clusters in the 
disk, 7 — 9 Gyr f3§ which are based on MSTO 
ages. This suggests that the age estimates based 
on the MSTO are reliable, and hence, that the 
GC age estimates are free of systematic errors. 

7. Summary 

Globular clusters are the oldest objects in the 
universe which can be dated. Absolute GC ages 
based on the luminosity of the MSTO are the 
most reliable JH) and lead to the tightest con- 
straints on the age of the universe. MSTO ages 
for GCs determined by a number of different re- 
searchers agree well with each, and have not ap- 
preciable changed for a number of years 0. A 
number of independent tests (summarized in 



of the stellar evolution models suggest that cur- 
rent stellar models are a good representation of 
actual stellar evolution, and hence, that there are 
no unknown systematic errors in the GC age esti- 
mates. A detailed study of the known uncertain- 
ties has lead to the conclusion that oldest GCs 
are 14. 6± 1.7 Gyr old, with a one-sided, 95% con- 
fidence limit that these clusters are older than 
12.1 Gyr To this age, one must add some es- 
timate for the time it took GCs to form after the 
Big Bang. Estimates for this formation time vary 
from 0.1-2 Gyr. To be conservative, the lower 
value is chosen. Thus, the universe must be older 
than 12.2 Gyr. This mimimum value of t requires 
that H D < 72km/s/Mpc if Q = 0.1, A = 0, or 
H Q < 54km/s/Mpc for a flat, dominated universe 
(CI = 1.0, A = 0). 

Acknowledgments 

Parts of this review are based on work done 
with P. Demarque, P. Kernan, and L. Krauss 
(Monte Carlo of GC ages), and Sarbani Basu and 
J. Christensen-Dalsgaard (helioseismology). I am 
grateful to Ata Sarajedini for furnishing Figure 1. 

REFERENCES 

1. Durrcll, P.R. & Harris, W.E. 1993, AJ, 105, 
1420 

2. Schwarzschild, M. 1958, Structure and Evo- 
lution of the Stars (Dover, New York) 

3. Kim, Y.-C, Fox,P.A., Sofia, S. & Demarque, 
P. 1995, ApJ, 442, 422 

4. For example, see: Langer, G.E., Kraft, R.P., 
Carbon, D.F., Friel, E. & Oke, J.B. 1986, 
PASP, 98, 473; and Kraft, R.P. 1994, PASP, 
106, 553 

5. The fact that the luminosity of the main se- 
quence turn-off is the most reliable age in- 
dicator for GCs is well established. See, for 
example Sandage, A. 1970, ApJ, 162, 841; 
Renzini, A. 1991, in Observational Tests of 
Cosmological Inflation, eds. T. Shanks, et ai, 
(Dordrecht: Kluwer), 131 

6. Carney, B.W., Storm, J. & Jones, R.V. 1992, 
ApJ, 386, 663 

7. For a discussion of the history of GC age es- 



10 



9. 

10 



11 



12 



13 



timates see: Demarque, P., Deliyannis, CP. 
& Sarajedini, A. 1991, in Observational Tests 
of Cosmological Inflation, eds. T. Shanks, et 
al. (Dordrecht: Kluwer), 111. For recent GC 
age estimates by various workers see: Berg- 
busch, P.A. & VandenBerg, D.A. 1992, ApJS, 
81, 163; Chaboyer, B., Sarajedini, A. & De- 
marque, P. 1992, ApJ, 394, 515; Salaris, M. 
Chieffi,A. & Straniero, O 1993, ApJ, 414 
580; Renzini, A. et al. 1996, |astro-ph/9604179| 
Chaboyer, B., Demarque, K., Kernan, P.J. & 
Krauss, L.M. 1996, Science, 271, 957 (CDKK) 
Chaboyer, B. 1995, ApJ, 444, L9 
Proffitt, C.R. & VandenBerg, D.A. 1991, 
ApJS, 77, 473; Chaboyer, B.C., Deliyannis, 
C. P., Demarque, P., Pinsonneault, M. H. & 
Sarajedini, A. 1992, ApJ, 388, 372 
Christensen-Dalsgaard, J., Proffitt, C.R. & 
Thompson, M.J. 1993, ApJ, 403, L75 
Chaboyer, B. & Demarque, P. 1994, ApJ, 433, 
510 

See, for example: Bahcall, J.N. 1989, Neu- 
trino Astrophysics (Cambridge: Cambridge 
U.P.); Bahcall, J.N. and Pinsonneault, M.H. 
1992, Rev.Mod.Phys. 64, 885 

14. For high temperature opacities see: Igle- 
sias, C.A. & Rogers, F.J. 1991, ApJ, 371, 
408; Seaton, M.J., Yan, Y., Mihalas, D. & 
Pradhan, A.K. 1994, MNRAS, 266, 805. For 
low temperature opacities see: Kurucz, R.L. 
1991, in Stellar Atmospheres: Beyond Classi- 
cal Models, ed. L. Crivellari, I. Hubeny, D.G. 
Hummer, (Dordrecht: Kluwer), 440; Neu- 
forge, C. 1994, in The Equation of State 
in Astrophysics, IAU Colloquium 147, eds. 
G. Chabrier and E. Schatzman (Cambridge: 
Cambridge U.P.), 618 

15. Rogers, F.J. 1994, in The Equation of State in 
Astrophysics, IAU Coll. 147, ed. G. Chabrier 
and E. Schatzman (Cambridge: Cambridge 
U.P.), 16 

16. Chaboyer, B. & Kim, Y.-C. 1995, ApJ, 454, 
767 

17. Mihalas, D., Hummer, D.G., Mihalas, B.W. 
& Dappen, W. 1990, ApJ, 350, 300 

18. Kernan, P.J. & Krauss, L.M. 1994, Phys. Rev. 
Lett., 72, 3309; Copi,C.J., Schramm, D.N. & 
Turner, M.S. 1995 Science, 267, 5195 



19. Lambert, D.L. 1989, in Cosmic Abundances 
of Matter, AIP conference proceedings 183, 
ed. C.J Waddington (New York: American 
Institute of Physics), 168 

20. See, for example: Bessell, M.S., Sutherland, 
R.S. & Ruan, K. 1991,ApJ, 383, L71; Tomkin, 
J., Lemke, M., Lambert, D.L. & Sneden, C. 
1992, AJ, 104, 1568; King, J.R. 1993, AJ, 106, 
1206; Nissen, P., Gustafsson, B., Edvardsson, 
B. & Gilmorc, G. 1994, A&A, 285, 440 

21. Walker, A.R. 1994, AJ, 108, 555 

22. Degrinnoccnti, S., Weiss, A. & Leone, L. 



1996, astro-ph/9602066 



23. Bergbusch, P.A. 1993, AJ, 106, 1024 

24. Faulkner, J. & Swenson, F.J. 1993, ApJ, 411, 
200 

25. Christensen-Dalsgaard, J. 1996, to appear in 
Proc. VI IAC Winter School: The Structure 
of the Sun, ed. T. Roca Cortes (Cambridge: 
Cambridge U.P.) 

26. For example, see: Basu, S. & Anita, H.M. 
1995, MNRAS, 276, 1402 

27. Basu, S., Antia, H.M. & Narasimha, D. 1994, 
MNRAS, 267, 209; Monteiro, M.J.P.F.G., 
Christensen-Dalsgaard, J. & Thompson, M.J. 
1994, A&A, 283, 247 

28. Bahcall, J.N. & Pinsonneault, M.H. 1995, 
Rev. Mod. Phys., 67, 781 

29. See, for example: Hata, N., Bludman, S., 
Langacker, P. 1994, Phys. Rev. D49, 3622 

30. Data from the recently completed Super- 
Kamiokanda, and soon to be completed SNO 
detector will provide a definitive test of 
the standard solar model. For example, see: 
Hata, N. 1994, to appear in Proceeding from 
the Solar Modeling Workshop, ed. B. Bal- 
antekin & J.N. Bahcall (World Scientific), 
hep-ph/9404329| 

31. Oswalt, T.D., Smith, J.A., Wood, M.A., 
Hintzen, P.M. & Sion, E.M. 1996, submited 
to Nature 

32. Garnavich, P.M. VandenBerg, D.A., Zurek, 
D.R. & Hesser, J.E. 1994, AJ, 107, 1097; 
Kaluzny, J. & Rucinski, S.M. 1995, A&AS, 
114, 1 



